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Abstract 

We calculate the contribution of the partonic processes gg — > WZqq and gg —* 
W^fqq to WZ and W^y pair production at hadron colliders, including anomalous triple 
gauge-boson couplings. We use the helicity method and include the decay of the W and 
Z-boson into leptons in the narrow-width approximation. In order to integrate over 
the qq final state phase space we use an extended version of the subtraction method to 
NNLO and remove collinear singularities explicitly. Due to the large gluon density at 
low x, the gluon induced terms of vector-boson pair production are expected to be the 
dominant NNLO QCD correction, relevant at LHC energies. However, we show that 
due to a cancellation they turn out to provide a rather small contribution, anticipating 
good stability for the perturbative expansion. 



1 Introduction 



A direct test of the non-abelian nature of the weak interactions can be made through 
the study of vector-boson pair production. These processes already involve triple gauge- 
boson couplings at tree level and, therefore, offer a promising framework for their study. 
In the Standard Model, the triple gauge-boson couplings are completely fixed by gauge 
symmetry. However, these couplings can be modified by new physics occurring at 
a higher energy scale. Given that the triple gauge-boson couplings have been much 
less precisely measured than, for example, the couplings of gauge bosons to fermions, 
it is natural to search for new physics by looking for anomalous triple gauge-boson 
couplings. 

Trilinear gauge-boson couplings have been studied at LEP2 [|IJ and the Tevatron 0. 
It should be emphasized that the studies at hadron colliders are complementary to 
studies at electron-positron colliders. While bounds on anomalous couplings obtained 
from LEP2 data are more stringent than those from hadron colliders, the latter offer 
the possibility to study these couplings at a higher centre of mass energy, where effects 
of anomalous couplings are enhanced. The amount of available data from hadron 
colliders will dramatically increase with the Run II at Tevatron and even more so with 
the LHC. This will allow the study of triple gauge-boson couplings to be taken one step 
further. Of course, all possible vector-boson pair production processes will be analysed. 
However, in this paper we concentrate on W ±r y and W ± Z pair production. 

In order to get reliable theoretical predictions for vector-boson pair production it 
is important to include higher-order corrections. One-loop QCD corrections, treating 
the vector bosons as stable particles, were computed for W ±r y production || and 
W^Z production M ten years ago. Since vector bosons are identified through their 
decay products it is essential to include the decay of the bosons into fermions. This 
allows arbitrary cuts to be added to the kinematics of the final states and therefore 
facilitates the comparison of theoretical predictions with experimental measurements. 
The calculation of these processes in the narrow-width approximation but retaining 
spin information via decay-angle correlations only in the real part of the amplitudes 
has been published in ||. Anomalous couplings for the W ±r y || and W ± Z processes 
have also been included. Finally, analytic results of the one-loop amplitudes for vector- 
boson pair production with the subsequent decay into lepton pairs in the narrow-width 
approximation have been presented in ||. In |ITJ| these amplitudes have been used 
to obtain numerical results including all decay-angle correlations. Steps beyond the 
inclusion of 0(a s ) corrections in the narrow-width approximation have been made 
PH , where non-doubly resonant diagrams have been included and in |12[], where 



in 



one-loop logarithmic electroweak corrections have been calculated. 

The size of the one-loop QCD corrections depends crucially on the physical quantity 
under investigation. Usually they are of the order of a few 10%, but in some particularly 
interesting cases, for example the high px tail of transverse momentum distributions 
they can increase the leading order cross section by several 100%. The reason for 
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such huge corrections is that the next-to-leading order correction to vector-boson pair 
production consists of two parts. First, there is the one- loop correction to the leading 
order partonic process qq — > VV. Second, there are new partonic channels, qg — > VVq, 
that contribute to the cross section at NLO. Even though these processes are suppressed 
with respect to the leading order partonic process qq — > VV by the strong coupling 
constant a s , they can easily be more important numerically, since they are enhanced by 
the gluonic parton distribution function. This is particularly the case for the LHC (or 
even a much higher energy Very Large Hadron Collider, or VLHC), where gluon induced 
processes become increasingly dominant. Unfortunately, the kinematical regimes where 
these contributions are largest are exactly those where effects from anomalous couplings 
are expected to manifest themselves. Therefore, it is important to gain a better control 
of the theoretical predictions in these kinematical regimes. 

One way to improve the situation is to impose a jet veto. This suppresses the effect 
of the new partonic channels that are opened at higher order [p], || [10[ . On the other 
hand, such a veto also reduces the amount of data considerably. Another possibility is 
to include even higher order terms. Of course, a full NNLO calculation of vector-boson 
pair production is not to be expected in the near future. Such a calculation would 
have to include two-loop corrections to the partonic subprocesses qq — > VV, one-loop 
corrections to qg — > VVq and, finally, (amongst others) processes with gg in the initial 
state. Motivated by the observation that loop corrections to partonic processes tend to 
be of the expected moderate size and huge corrections are mainly due to the opening of 
new channels, in this paper we include all NNLO terms that are maximally enhanced 
by the gluon distribution functions. 

In general, there are two classes of such contributions. First, there is the process 
gg — > VV. These processes have been studied previously [pT3|| . The amplitude for 
this process vanishes at tree level but may be non-zero at one-loop. In the case of 
W ±r y and W^Z pair production, however, the amplitude for this process vanishes at 
all orders, due to charge conservation. Second, there are the processes gg — > VVqq. 
The amplitudes for these processes have been calculated in |14| and have been used 
to compute cross sections for the production of a vector boson pair together with 
two jets. We extend this work by not requiring any jets to be observed in the final 
state. Hence, the amplitudes gg — > VVqq have to be integrated over the qq final state 
phase space. This results in infrared singularities that will have to be absorbed in 
the parton distribution functions. Furthermore, we also extend the previous work by 
adding anomalous couplings. 

We find that the contribution of the gluon induced processes to WZ and Wj 
production is smaller than may have been anticipated, due to a sign change in the 
hard scattering part. We also find that this NNLO term is affected less by anomalous 
couplings than the corresponding LO and NLO terms. 
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2 Calculation 



Vector boson pair production has been calculated up to NLO in the helicity method. 
However, in an NLO calculation for pp — > WZ, whilst the process qq — > WZ is included 
at NLO, the process qg — » VTZg is only calculated at LO. This means that we only have 
a leading order calculation in the circumstances where gluon-induced subprocesses are 
important, as they are expected to be at the LHC due to the large gluon density at 
low x. Therefore we want to include the next order in a s , the process gg — > WZqq 
(and similarly for Wj). 

For W~Z production we will use the labelling g 1 g 2 13 v± 1' 5 1' 6 q 7 q s , where g\ and g 2 
are the incoming gluons, leptons 3 and 4 are the decay products of the W and 5 and 6 
are from the Z. Clearly, Wj will have one less particle but this is a trivial modification. 
It is also straightforward to transform the calculation for other vector bosons. 

In calculating gg — > WZ, we use a version of the subtraction method |T3| to cancel 



the infrared singularities analytically and evaluate the finite remainders numerically. 
We show the finite parts of the cross section in equations (^j), ( |ITD and flT5|). 
A generic differential cross section is written as: 



da AB (p A ,p B ) =S/ dx i dx 2 fa/A(x A )f b/B (x B ) da ab (x A p A ,x B p B ) 

ab 



where A and B are hadrons while a and b are partons and f a / A and f / B are the 
parton distribution functions. The subtracted partonic cross section, da, is finite for 
any infrared safe observable. 

In a general case, in order to obtain da, soft and final state collinear singulari- 
ties have to be cancelled with the singularities coming from virtual corrections. The 
remaining initial state collinear singularities have to be factored out into the parton 
distribution functions. In an extension of the subtraction method to NNLO, this de- 
mands the subtraction of all possible singularities as well as their analytical integration 
in the corresponding limits. Up to now, such a general framework has not been devel- 
oped. Only a limited number of analytical calculations to NNLO accuracy have been 
performed. 

The main difficulty for the implementation of a subtraction procedure at NNLO is 
the appearance of several new (multiple parton) kinematical configurations where soft 
and/or collinear singularities arise. These are: the triple collinear case (three partons 
become collinear), the double soft case (two partons are soft), the soft-collinear case 
(one parton is soft and the other collinear) and the double collinear case (two pairs of 
partons are independently collinear). The single soft and collinear singularities, which 
appear also at NLO, have to be added to that list. 

Even though the behaviour of tree and one-loop amplitudes in those limits has been 



recently obtained |L6[ , it is not yet clear how to unify all singularities in a single subtrac- 
tion term, or how to subtract all of them independently while avoiding double counting 
and allowing the analytical integration of the subtracted terms to be performed. 
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Our case, however, is particularly simple. Since we include only gg induced pro- 
cesses at NNLO there are no final state collinear singularities (see below about final 
state singularities in the W'j case). There are also no soft singularities, since there are 
no gluons (or qq pairs) in the final state. Therefore, the only singularities that arise 
are the single collinear ones, which can be treated in a similar way as a typical NLO 
calculation, and the double collinear ones, i.e., the collinear splitting of each initial 
state gluon into a qq pair, a 'genuine' NNLO contribution. 

In order to absorb the remaining initial state collinear singularities we use 

fa/d = 6*6(1 -x) - ~P ad (x,0) - f^) 2 Q ad ( x ) . (2) 

This is the NNLO expression in MS and as usual we define 4 = - — 7e + log47r. 
P a d{x, 0) is the one-loop Altarelli-Parisi splitting function for e = (4 dimensions). 
The 0(a 2 s ) term Q ad , including contributions from the two-loop splitting functions, 
does not actually contribute in our case. 

Comparing subtracted and unsubtracted terms and using eq. (0), we obtain the 
relation: 



dafj(k u k 2 )=da^(k 1 ,k 2 ) 
+ ^ J dxi 0zP qg {xi) \ da^(xik u k 2 ) + ^ J dx x 0=Pq g (xi) \ daj$ \x x k u k 2 
J dx 2 {^zP qg {x 2 )\ da^' '(&!, x 2 k 2 ) + ^ J dx 2 {^zPqg{x 2 ) \ da^ (ki, x 2 k 2 



2tt 



+ 

i 

+ (^) J dxidx 2 (^P qg (xi)^ (^Pq g (x 2 )^ daf^(xiki,x 2 k 2 ) 

+ (^) J dxidx 2 CzPqgixx^j 0zP qg {x 2 )\ daf q \xxk Xl x 2 k 2 ) , (3) 

where fci, k 2 are the momenta of the incoming gluons 1 and 2. The superscript (2) 
denotes the NNLO term, while (1) is NLO and (0) is leading order. In eq. (^) all terms 
on the right hand side are separately divergent but the sum is finite. 

Once we have this expression, we can write out the counterterms in terms of en- 



ergy and angle variables and cancel the poles explicitly. Following |T5| we use the 
parameterisation 



h = ^(1,0,1) (4) 



h = ^(1,0,-1) (5) 



h = & ^1, V 1 -Vi&r,Vi) , ( 6 ) 

where ki is the momentum of an outgoing (anti) quark, ct is a unit vector in transverse 
momentum space, yi = cos#j is the angle variable — 1 < yi < 1 and £j is the rescaled 
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energy variable < £j < 1. Performing part of the phase-space integration analytically 
we can rewrite eq. ([3]) as a sum of three separately finite pieces 



(7) 



with terms da^ n ' 8 ' 7 ' 6 ^ denoting the 8,7 and 6-parton finite parts. The explicit form of 

^,j(fi n ,8) j g giy en 



da {jin ' : 



;i - y 2 7 )(l - jfi)M<$ (h, k 2 , {h} 3;6 , k 7 , k 8 ) -j^frfcp^) V(y s) 



64(2 



d£ 7 d^s dy-j dy 8 dip 7 d<p 8 d$ 



(3-6) , 



(8) 



where d$( 3 _ 6 ) denotes the phase-space integration over the vector-boson decay prod- 
ucts. M.gg denotes the squared amplitude summed (averaged) over helicities, including 
the flux factor l/(2si 2 ). V{yi) is given by: 



^w- ^[(1^)^(1^)^ 



(9) 



The distributions (j^jp) have been introduced in [T5\ and are defined for an arbitrary 
test function f(yi) through 



J(Vi)) 



1 d f(yi)-f(n)o&yi-i + 6i) 

-1 1 ± Vi 



(10) 



The two V distributions in eq.(|8]) perform the subtraction of both double and single 
collinear singularities of NNLO amplitudes, leaving a (numerically) integrable remnant. 
The 7-parton finite part reads 



{m£ ] ((1 - k 2 , {ki} 3fi , k 8 ) + M { p q (h, (1 - &)k 2 , {ki} 3fi , k 8 ) } 



a 



fs V{y%) d& d£ s dy 8 dip 8 d$ 



(3-6) 



+ g (£ 8 P<(1 - &) - P'<{1 - fc)) (1 - y 2 7 ) 
{M$ ((1 - k 2 , {ki} 3fii k 7 ) + M%j (h, (1 - Cs)^, {kh, 6 , k 7 )} 



S12 



8(2tt) 3 
where we introduced 



^V{y 7 ) d£ 7 d£ 8 dy 7 dip 7 d<$> 



(3-6) , 



'11} 



C 7 = loj 



gig Si £? 
2/i 2 



Cr = lot 



S12 Si £g 
2/j, 2 



(12) 



and split up the unregularized Altarelli-Parisi function into the 4- dimensional piece, 
P^d, and the piece proportional to e, P^. For completeness we give the explicit form 

^£(1-6) = i(l-26 + 2e?) (13) 

P£{i-&) = 6(6-1)- (14) 

In this case, the single collinear singularities of the NLO amplitudes are subtracted by 
the appearance of a single V distribution. 

Finally, the 6-parton finite piece is given by 

= (£) 2 (ap<(i - w - ^(i - (^(i - &) - n?a - &)) 

{A^ ((1 - &)M1 " £8)M^}(3,6)) 

+ ((1 - &)&!, (1 - £ 7 )A; 2 , M {3 ,6)) } d6 ^(3-6) • (15) 

We note that 5j is an arbitrary quantity and even though da^ n ' 8 \da^ n '^ and 
rfcr( fin > 6 ) depend on it, this dependence cancels exactly in da^J. Hence, making sure 
that the numerical results are independent of 5i is a useful check for our Monte Carlo 
program. 

In calculating the cross section as above, we use 6 parton and 7 parton amplitudes 
as given in ||. We also require the 8 and 7 parton amplitudes for gg — > qqWZ and 
99 ~ ¥ qqW'j respectively. These amplitudes have been calculated previously fT3fl . We 
recalculated these amplitudes using the helicity method and including anomalous cou- 
plings. For the triple gauge vertex for W£ {p+)W~{p-)pV tl (q), including the anomalous 
terms, we use 



1 1 „v i .y i \v P 



+ ( q -p_)f>g«»l( J gV+ K V+\ 



2 



Ml' 



1 „ A" X v 



^-P-n-9-''- 2 f M + W ^ q '% (16) 

where V can be Z or 7 and all momenta are taken to be outgoing. In the Standard 
Model, g\ = k v = 1 and = 0. In order to preserve electromagnetic gauge invariance 
we always set gj to the Standard Model value 1. 

We also make use of form factors to avoid the violation of unitarity that would 
otherwise result from this vertex at high energies. The form factors used is are the 
conventional ones: 

v AgY v Ak v v AX V 

Adl ~* (1 + s/A 2 ) 2 ' ~* (1 + S/A 2 ) 2 ' (1 + s/A 2 ) 2 { } 
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where A is the scale of the new physics that causes the anomalous couplings. 

Photons with large transverse momentum can be produced in hadronic collisions not 
only directly, but also, and significantly, from the fragmentation of a final state parton. 
A full perturbative calculation of a process involving the production of photons should 
in principle include the calculation of both direct and fragmentation components, since 
only their sum is physically well defined beyond LO. A calculation of the fragmentation 
part is not even available at NLO for the process of interest in this paper. Fortunately 
there is a way to suppress the fragmentation contribution, which actually constitutes 
a background to the search of anomalous couplings, by requiring the photons to be 
isolated from the hadrons. The usual way to perform the isolation is to require the 
transverse hadronic momentum in a cone around the photon to be smaller than a 
fraction of the transverse momentum of the photon. In this way, the contribution 
of the fragmentation component can be reduced to the percent level. In this paper, 
we will use the isolation procedure introduced by Frixione |l7j which allows us to 
completely suppress the fragmentation component. Therefore, we reject all events 
unless the transverse hadronic momentum deposited in a cone of size i?o around the 
momentum of the photon fulfills the following condition 

Y,PTi9(R - R iy ) < PTl ( }~° 0S * ) , (18) 
i VI - cosR J 

for all R < Rq, where the 'distance' in pseudorapidity and azimuthal angle is defined 
by i?j 7 = \J (r/i — r/ 7 ) 2 + (4>i — </> 7 ) 2 . In this way, only soft partons can be emitted 
collinearly to the photon in the direct contribution and, therefore, no final state quark- 
photon collinear singularity arises. For the purposes of our computation this allows us 
to perform the calculation of the gg initiated contribution at NNLO without needing 
to do any subtraction of the corresponding singularity. While the possibility of per- 
forming such isolation at the experimental level is under investigation, the choice of 
this particular method will not alter at all the conclusions of our calculation in the W'y 
channel. Particularly, we have checked that the results at NLO using the procedure 
in eq.([TJ) are very close to the ones obtained performing the usual 'cone' isolation 
procedure fLOf . 

In this calculation, contributions from b and t quarks have been neglected. It is 
assumed that these will be suppressed by the large top quark mass. This might not be 
a particularly good approximation for large energies, but in view of the smallness of 
the gluon induced corrections a more detailed treatment doesn't seem to be justified. 



3 Numerical Results 

For the numerical results presented in this section we use the MRST 2001 parton 
distribution functions [TTU] with the one-loop expression for the coupling constant 



(a s (Mz) = 0.119). The factorization and renormalization scales are fixed to hf = vs, 
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fl>R 



+ \ {prw +Ptz) i n the case of WZ production and to fip 



fin = yMjy + Pt-j m the W'y case. Notice that in order to compute the gg induced 
contribution we also use the same NLO combination of parton distribution function 
and coupling constant. The effect of changing to NNLO parton distributions, not fully 
available yet, is expected to be small and will not alter the conclusion of our analysis. 

The masses of the vector bosons have been set to Mz = 91.187 GeV and Mw = 
80.41 GeV. We do not include any electroweak corrections, but choose the coupling 
constants a and sin 2 6y/ in the spirit of the improved Born approximation [2D||. For the 
couplings of the vector bosons with the quarks we use a = a(M z ) = 1/128 whereas 
for the photon coupling we use a = 1/137. We neglect contributions from initial state 
top quarks and use the following values for the Cabibbo-Kobayashi-Maskawa matrix 
elements: \V u d\ = \V CS \ = 0.975 and |V^ S | = \V c d\ = 0.222. Note that we do not include 
the branching ratios for the decay of the vector bosons into leptons. 

For all plots we use a set of standard cuts. For charged leptons we require px > 
20 GeV and 77 < 2.5. In addition, we require a missing transverse momentum p™ lss > 
20 GeV. The photon transverse momentum cut we use is > 20 GeV, while for the 
isolation prescription in eq. (113) we set Rq = 1. 
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Figure 1: Partonic cross sections for W'y: a q q, a qg and a, 



In Figures [I] and |2] we show the MS subtracted partonic cross sections for the 
different initial states qq, qg and gg. These are produced by binning the results at 14 
TeV by y/I. It can be seen that the outcomes for W'y and WZ are very similar. For 
better visibility a gg has been increased by a factor 10 in the plot. 

As can be observed, there is a reasonably good perturbative convergence at the 
level of the partonic cross section, considering that a q q is of 0(1 + a s ), a qg of 0(a s ), 
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Figure 2: Partonic cross sections for WZ : a q q, a qg and a gg 

and <jgg of 0(a 2 s ), with a s ~ 0.1 for the typical scales of these processes. In the MS 
both the qq and qg contributions are positive in the whole range of s, while the new 
piece, the gg partonic cross section becomes negative only for small values of s, close 
to the threshold for the production of the vector bosons. 

The situation changes considerably when the physical hadronic cross section is 
computed, shown in Figures [3] and |j. In Figure [3] we plot the transverse momentum 
distribution of the photon, whereas in Figure £| we see the pt of the lepton produced 
by the decay of the W. Again the form of the distributions is very similar for the WZ 
as compared to the case. 

As indicated in the introduction, because of the large qg luminosity, the hadronic 
contribution due to this initial state becomes even larger than the 'leading' qq contri- 
bution. This is particularly clear at large transverse momentum, where the relative 
increase in the luminosity exceeds the effect of the suppression due to the extra a s 
coupling observed at the partonic level. This is not the case for the gg contribution, 
which remains rather small, and for the observable studied here provides a negative 
contribution to the hadronic cross section. 

The features of both qg and gg contributions to WZ and W'y production can be 
qualitatively understood in the following way: the qg and gg luminosities are steep 
functions of the momentum fraction carried by the partons, mostly due to the fast 
increase of the gluon distribution g(x) when x decreases. Even though the hadronic 
center of mass energy S is very large, the average value of the partonic one (s) = 
(xiX2)S can be considerably smaller, and actually closer to the minimum needed to 
produce the gauge bosons with the required transverse momentum. Therefore in the 
convolution between the parton distributions and the partonic cross sections shown in 
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Figure 3: W'y production: p T distribution for LHC, qq, qg and gg pieces separately 

Figures and 0, the hadronic result will mostly pick up the features of the partonic 
cross section at low values of s. In the qg case, the partonic cross section is positive 
and has its maximum at low s, giving as a result a large hadronic contribution. In 
the gg case, the partonic cross section has a change of sign at low values of I, with 
the corresponding compensation between negative and positive contributions when 
convoluting with the parton densities. Therefore the gg hadronic contribution results 
in a negative and small correction to the NLO cross section. 
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Figure 4: WZ production: px distribution for LHC, qq, qg and gg pieces separately 
In order to investigate whether the inclusion of anomalous couplings changes this 
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picture substantially, we show in Figure |5| the partonic cross sections with the following 
anomalous couplings: 

g[ = 1; gf = 1.13; K i = 1.2; k z = 1.07; A 7 = \ z = 0.1 (19) 

We use the form factors as given in eq. ( |ITD and we take A, the scale of new physics, to 
be 2 TeV. Figure |5] is given with the same scale as Figure [| for comparison. It is clear 
that, with the values for the couplings that we have chosen, we do not see a significant 
enhancement of the gluon induced term. In fact, the gluon induced corrections become 
even less relevant, as the qg and especially the qq terms are increased substantially. 



20 


>- da xl0 6 [pb/GoV]\ 








; dVs" \ 


solid qq 








dashed qg 




15 




\ dotted gg (xlO) 




10 








5 













i , , , , i , , , , 





2500 5000 7500 10000 12500 

Vs [GeV] 



Figure 5: Partonic cross section a qq -, a qg and a gg for WZ production with anomalous 
couplings as given in eq. fllQ) 

In addition to these results, we also considered a hypothetical VLHC or Very Large 
Hadron Collider, which would run at a centre of mass energy of 200 TeV, to see whether 
the gg part would become significant at the higher energy (cuts etc. all remained the 
same as previously). It was found that the gg part did not become more important. 
The contributions from the gg term remained at the 1% level. 

4 Conclusions 

We have presented expressions for the finite cross section for vector-boson pair produc- 
tion from a gluon-gluon initial state. We have then implemented these terms and the 
appropriate matrix elements into a Monte Carlo integration. 

We find the contribution of the gluon induced terms in WZ and W'j production 
to be surprisingly small. Even the addition of anomalous couplings enhances the gg 
induced term substantially less than the qq and qg terms. We have shown that this is 
due to a change of sign in the hard scattering. 
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This is a good result from the point of view of experimental predictions of vector- 
boson pair production, at least in the WZ and cases, as the impact of NNLO 
terms does not appear to be particularly significant. It also seems that the gg term 
will not be amplified excessively by anomalous couplings, though we have only taken 
an example of anomalous couplings and have not made a detailed study. 
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